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Abstract 


The  dependence  of  the  stress  field  on  a  continuous  periodically  varying 


stiffness  which  approximates  a  layered  material  is  studied.  This  dependence 

a 

is  compared  to  the  solution  for  the  usual  approximation  of  'converting  this 
class  of  problems  to  the  superposition  of  an  anisotropic  one  and  a  residual 
boundary  problem.  The  region  of  applicability  of  the  usual  approximation  is 
found  to  be  as  expected  geometrically  but  surprisingly  large  as  far  as  variation 


of  stiffness 


is  concerned . 
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Introduction 


Inhomogeneities  in  stiffness  of  isotropic  elastic  materials  occur  in 
several  areas  of  technology.  In  particular  various  applications  in  rock  mechan¬ 
ics,  layered  materials  and  composite  materials  are  related  directly  to  such 
inhomogeneities  [1,2,3]*.  When  an  elastic  modulus  varies  periodically  with  one 
Cartesian  coordinate  the  usual  situations  occuring  in  rock  mechanics  and  layered 
materials  can  be  represented.  This  work  is  concerned  with  such  a  case. 

Let  the  wavelength  of  elastic  modulus  variation  be  L  and  some  charac¬ 
teristic  dimension  for  a  particular  boundary  value  problem  be  c  .  When 
L  «  c  an  analytical  procedure  has  been  described  [2]  for  "converting1'  the 
inhomogenous  problem  into  the  superposition  of  two  problems  which  are  1)  an 
anisotropic  material  problem  with  the  original  configuration  and  boundary  condi¬ 
tions  and  2)  a  residual  boundary  problem  whose  effect  on  the  stress  distribution 
should  be  important  only  a  few  wavelengths  into  the  body  from  the  surface. 

This  superposition  technique  is  currently  the  standard  procedure  for 
dealing  with  this  class  of  problems.  The  reason  is  that  more  solutions  and  tech¬ 
niques  [4]  are  available  for  solving  anisotropic  problems  than  for  inhomogenous 
problems.  Also,  the  residual  boundary  problem  can  be  approximated  as  a  local 
one  independent  of  the  configuration  of  the  particular  problem  under  considera¬ 
tion.  Presumably  this  inhomogenous  residual  boundary  problem  can  be  solved 
"once  and  for  all"  for  a  particular  type  of  stiffness  variation  and  inclination 
of  the  boundary  to  the  layers. 

The  residual  boundary  problem  has  received  some  attention  in  the  liter¬ 
ature.  When  the  variation  of  the  stiffness  is  discontinuous  (e.g.,  bonded, 
layered  materials)  the  residual  boundary  problem  has  been  considered  by 

*Numbers  in  brackets  refer  to  the  References  at  the  end  of  the  paper. 
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Bogy  [5],  When  the  stiffness  variation  is  sinusoidal  with  small  amplitude  com¬ 
pared  to  the  mean  and  the  boundary  is  perpendicular  to  the  layers  the  solution 
to  the  residual  boundary  problem  may  be  found  in  [6].  In  [5,6]  a  strong  depend¬ 
ence  of  the  stresses  on  Poisson's  ratio  is  found. 

The  object  of  this  work  was  to  solve  an  elementary  problem  with  contin¬ 
uously  varying,  largo ,  stiffness  changes  without  resorting  to  the  superposition 
technique  described  above.  By  solving  this  inhomogenous  elastic  problem  of  an 
infinite  slab  tinder  tension  perpendicular  to  the  layers  it  was  hoped  the  region 
of  applicability  of  the  superposition  technique  could  be  more  firmly  established. 
It  is  best  to  have  an  exact  solution  to  an  appropriate  inhomogenous  problem  for 
this  purpose  but  none  are  currently  available. 

Figure  1  shows  the  problem  considered  here.  The  infinite  slab  has  a 
thickness  2c  and  the  wavelength  for  the  shear  modulus  variation  is  L  .  The 
problem  is  two  dimensional  so  that  a  very  elementary  numerical  procedure  may  be 
used  to  approximate  values  for  the  stresses.  The  solution  for  the  dimensionless 
stresses  ,  a^/o^  and  0^/0  AVG  depends  on  the  following  constant 

parameters : 

v  *  Poisson's  ratio 

L/c  —  Material  property  variation  wavelength /characteristic  dimension  of  the  body 
e  =  A  measure  of  the  stiffness  variation 
The  practical  ranges  of  these  parameters  chosen  here  for  study  are 
0.1  {  v  {  0.5 
0.1  $  L/c  $  3.0 
0  (  c  $  0.7 

In  the  following  two  sections  the  problem  is  formulated  and  an  approxi¬ 
mate  solution  technique  outlined.  Then  the  solution  presented  in  [6]  which 
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employs  the  superposition  of  anisotropic  and  residual  boundary  solutions  is 
given  for  the  problem  of  Figure  1.  Subsequently  numerical  results  are  presented 
comparing  the  approximate  solution  obtained  here  to  the  perturbation  solution 
from  [6]  and  conclusions  are  drawn. 


Formulation  of  the  Problem 

Since  the  shear  modulus  G(x)  is  periodic  in  x  the  stresses  are 

also  periodic.  With  the  notation 


f(x)  =  (1  -  ecos(2irx/L)}  ^ 


(1) 


the  relevant  plane  strain  stress-strain  relations  are 
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where  u  and  v  are  the  x  and  y  displacements.  The  stresses 

and  a  can  be  deduced  from  the  Airy  stress  function  ij»(x,y)  so  that  they 
xy 

satisfy  the  force  equilibrium  equations  for  this  static  plane  strain  problem  as 


xx 


yy 


yy 


XX 


and 


xy 


=  *, 


xy 


(3) 


while  o  is  given  by  v(o  +  o  )  and  o  as  well  as  o  vanish, 
zz  6  J  xx  yy  xz  yz 

Combining  equations  (2)  and  (3)  to  eliminate  u  and  v  yields 


f’x 

g('l')  -  ^»xxxx  +  2^’Xxyy  +  ^*yyyy  +  2  f  ^’xx  +  ^’yy^’x 


j. -  (ip,  -  T-^—  l//,  )  =  0 

f  xx  1-v  yy 


(4) 
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Appropriate  boundary  conditions  for  this  problem  are  deduced  as  follows: 


°*y  *  °yy  vanish  m  y  s  ±c  *  'f>,xK‘x,±c)  ~  0  *  ^xy***10*  =  0 

no  net  shear  force  on  cross-section  ■*  4i>x(x,-c)  =  \|»,x(x,+c) 

(5) 

no  net  moment  on  cross-section  -*■  c{^,y(x,c)  +  ^,y(x,-c)}  =  i|»(x,c)-t|»(x,-c) 
net  axial  force/depth  is  2oAVQc  -*•  ^,y(x,c)  -  <l»,y(x,-c)  =  2oAV(,c 

It  is  noted  that  the  usual  step  of  replacing  v  in  the  plane  strain  formulation 
by  will  yield  the  formulation  for  plane  stress. 


Outline  of  Solution 

The  following  series  approximation  for  f  in  terms  of  the  (m-l)n 

arbitrary  constants  H.  .  (i  =  2,3... m;  j  =  l,2...n)  was  chosen 

J 


°AVfi 


(v2 


m/2 


I  H  .  (cos(irky/c)  -  (-l)k}  cos(2irqx/L) 
q=l  Zk*q 


m/2-1  n  . 

+  11  H_.  .  {cos[(2k+l)iry/2c]  -  (-1)  (2k+l)cos(iry/2c)}  cos(2irqx/L) 

k=l  q=l  ZK+'L;q 
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(6) 


from  equation  (6)  satisfies  all  the  boundary  conditions  given  in  equation  (5). 

The  H.  are  determined  using  the  h.  (x,y)  as  weighting  functions  so  that 
*.;q  x.;q 

L  +c 


h£.q  g(<l»)  dx  dy  =  0 


£  =  2,3. ..m;  q  =  1,2. ..n 


(7) 


x=0  y=-c 

when  ip  from  equation  (6)  is  substituted  into  (7).  The  resulting  (m-l)n 

2 

simultaneous  equations  for  tb«  H.../L  oAVQ  depend  on  L  ,  e  and  v  and,  of 
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course ,  m  and  n  .  These  equations  were  programmed  and  solved  on  a  digital 


computer  for  about  150  specific  choices  of  L,  e,  v,  m  and  n 


The 


A 

H._  ,/L  were  subsequently  used  to  find  °£j/°AyG  by  means  of  equations  (3) 

and  (6).  In  addition  to  the  dimensionless  stresses  °xx/°Av g’  °yy^0AVG  aTK* 

0  /o.„_  the  dimensionless,  plane  strain,  octahedral,  shear  stress  t/o  was 

xy  AVG  AVG 

evaluated  from 

T'«AVG  =  ^  ^°yy/0AVG^^  - 


+  <VW*>  1/2 


(8) 


These  dimensionless  stresses  were  evaluated  for  combinations  of  x/L  and  y/c 
of  interest.  Values  of  m  and  n  were  increased  in  trial  runs  and  correspond¬ 
ing  variations  in  t  /o for  fixed  L,  e  and  v  were  observed.  It  appears 
that  an  accuracy  of  a  few  percent  can  be  expected  using  m  and  n  equal  to  4 
in  the  range  of  L,  e  and  v  considered.  Due  to  the  large  number  of  cases 
computed  it  was  necessary  to  use  the  minimum  values  of  m  and  n  which  yield 
results  of  acceptable  accuracy. 


Perturbation  Solution 

Nusayr  and  Paslay  [6]  recently  determined  the  stresses  near  the  boundary 
for  the  material  described  in  Figure  1  using  a  perturbation  technique.  The 
stresses  near  the  boundary  in  the  present  problem,  therefore,  should  match  the 
solution  in  [6]  when  the  restrictions  of  [6]  apply.  These  restrictions  are  that 
L/c  and  e  be  small  compared  to  one.  The  solution  is 
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o  /o  =  l-e{v/(l-v)Hl-(2irc/L)[l-(y/c)]}e"(2irc/L)Cl"<y/c)]cos(2i»x/L) 

XX  AVu 

a  /ot_  =  -c{v/(l-v)Hl+(2irc/L)[l-(y/c)3}e"(2irc/L)Cl"(y/c):icos(2irx/L)  (9) 

yy  Avb 

°xy/oAVG  =  +e<v/(l.v)>(2xc/L)ri.(y/c)]e-(2-/L)Cl-(y/c)]sin(2wx/L) 

where  e  «  1  ,  L/c  «  1  ,  0  j  x  {  L/2  and  0(yjc.  In  this  paper 

equations  (9)  are  referred  to  as  the  perturbation  solution. 

Numerical  Results 

In  the  above  it  was  pointed  out  that  attention  here  was  restricted  to 
ranges  of  the  parameters  of  0.1$  0.5,  0  <  £  {  0.7  and  0.1  $  L/c  £  3.0. 

These  ranges  cover  the  region  of  most  practical  interest.  Computations  were 
made  for  enough  combinations  of  these  parameters  to  indicate  trends  in  the 
solutions  and  to  provide  curves  suitable  for  applications. 

Frequently  the  most  important  quantity  for  a  designer  to  have  is  the 
maximum  value  of  a  yield  criterion.  In  this  case  the  plane  strain  octahedral 
shear  stress  criterion  was  evaluated  and  designated  t  .  For  each  computed 
solution  the  maximum  value  of  T/aAVG  was  chosen  from  a  gr^d  of  comPuted  values 
over  the  region  0  *  x  $  L/2  and  0  $  y  $  c  .  In  every  case  the  maximum  value 
of  t/o4  _  was  on  the  x  =  L/2  boundary  of  this  region.  Figure  2  shows 
results  obtained  from  the  computations.  The  three  plots  for  different  values 
of  v  show  the  strong  dependence,  mentioned  above,  of  the  solution  on  Poisson’s 
ratio.  When  v  =  0.1  the  influence  of  the  considered  e  and  L/c  variations 
on  t  /o.„„  is  negligible.  When  v  =  0.5  the  influence  of  increasing  e 

from  0  to  0.7  can  almost  double  the  value  for  the  yield  criterion  Tmax/°AvG* 
The  influence  of  increasing  L/c  is  to  decrease  the  severity  of  the  change  of 


t  /o..,_  with  c  as  one  intuitively  might  expect.  The  curve  for  v  =  0.3 
max  AV6 

shows  behavior  intermediate  to  those  for  0.1  and  0.5  .  In  Figure  2  the 

i 

dashed  curves  show  the  predictions  of  the  perturbation  solution,  equations  (9). 

As  long  as  L/c  $  1.0  the  perturbation  solution  gives  art  acceptably  accurate  > 
prediction  for  t^/o^  .  ' 

An  effect  in  this  solution  that  can  be  easily  comprehended  is  the  devel- 

i  .  i 

opment  of  a  due  to  differential  lateral  contraction  at  x  =  0  and  x  =  L/2  . 

yy 

l 

A  tensile  o^/o^  develops  at  x  =  0  ,  y  =  0  while  a  compressive  o^/o^ 

developes  at  x  =  L/2  ,  y  =  0  .  These  are  the  minimum  and  maximum  values  of 

o  /o...„  in  the  region  0  $  x  $  L/2  and  0  $  y  S  c  .  In  each  computed  case 
yy  AVG 

the  value  of  the  maximum  o  was  greated  than  the  absolute  value  of  the  minimum 

yy 


.  The  perturbation  solution  gives  max 


o  .  .  Figure  3 . shows 

1  yy  min'  6 


results  for  a  in  the  same  format  as  Figure  2.  It  is  apparent  from  these 

yy  max 

curves  that  the  Poisson's  ratio  influence  is  present,  that  mgx  is  the  order 

of  2evo.„_  for  L/c  <  1.0  and  that  the  perturbation  solution  gives  reasonable 
AVG 

estimates  for  o  over  the  range  of  parameters  considered.  ‘ 

yy  max 


Conclusions  , 

.  i 

The  results  of  the  analysis  presented  here  help  to  make  clear  the  influ¬ 
ence  of  isotropic  inhomogeneities  in  isotropic  elastic  materials,  the  numerical 
results  in  Figures  2  and  3  indicate  the  influence  on  t^/o^  and 

a  /a  _  of  Poisson's  ratio  v  ,  the  inhcmogeneity  wavelength-characteristic 

yy  max  AVG 

length  ratio  L/c  and  the  inhomogeneity  strength  parameter  e  .  The  applicabil¬ 
ity  of  a  perturbation  solution  which  approaches  the  solution  to  the  problem  as 
L/c  -*•  0  and  e  -»•  0  is  shown  to  have  a  wider  applicability  in  the  solutiort 
to  the  problem  considered  here  than  the  author  originally  anticipated.  It  is  not 

l  '  ■ 

I 
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surprising  that  the  perturbation  solution  fails  to  adequately  predict  t  /o. vr 

max  AVo 

i 

for  L/c  >  1  however  the  range  of  e  values  for  which  t  /o.„_  and 

!  max  AVb 

Oyy  max/°^vG  are»  ^or  most  purposes ,  adequately  predicted  is  surprisingly 
large.  1 

!  S 
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Figure  Captions 


Figure  1 


Figure  2 


Figure  3 


Configuration  and  description  of  material  properties  for  the 
problem  under  consideration. 

Dependence  of  the  maximum  plane  strain  octahedral  shear  stress 

t  /o.„_  on  the  material  variation  parameter  e  ,  the  material 
max  AVG  r 

wavelength /characteristic  length  parameter  L/c  and  Poisson's 
ratio  v  .  Solid  lines  give  the  solution  presented  here  while 
dashed  lines  give  the  perturbation  solution  from  [6]  which  is 
independent  of  L/c  . 

Dependence  of  the  maximum  la  I  ,  o  .  on  the  material 

yy  yy  max 

variation  parameter  c  ,  the  material  wavelength/characteristic 
length  parameter  L/c  and  Poisson's  ratio  v  .  Solid  lines  give 
the  solution  presented  here  while  dashed  lines  give  the  perturba¬ 
tion  solution  from  [6]. 
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G(x)  =  SHEAR  MODULUS  =  G0  (l-€  cos  Zwj-)t  0<€<  I 
v  =  CONSTANT  POISSON'S  RATIO 
crAV6  =  AVERAGE  cr„  STRESS  ON  CROSS-SECTION 
PLANE  STRAIN 
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